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Inclusive production oiTJ/ip pair in proton-proton interation at LHCb 
is considered. This process is forbidden at leading order of perturbation 
theory, so such channels as double parton scattering, XbXc pair production 
with subsequent radiative decays of P-wave quarkonia, contributions of 
color-octet states, and NLO corrections are studied in details. For all 
these channels we present theoretical predictions of total cross sections at 
LHCb and distributions over different kinematical variables. According 
to presented in the paper results, double parton interaction gives main 
contribution to the cross section of the considered reaction. 
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1 Introduction 


Inclusive double production of heavy quarkonia was observed in NAS experiment, 
where production of two J/ijj mesons in ttN interaction was studied. The cross section 
of double J/'ip production was suppressed by about three orders of magnitude in 
comparison with single J/^j. Such suppression can be easily explained at leading order 
(LO) of perturbative QCD [H [2]- In low energy hadronic collisions quark-antiquark 
subprocess qq —)■ J/if) J/'ijj gives main contribution to the considered reaction, while in 
the case of high energy interactions gluonic subprocesses gg —?• are dominant. 

Leading order O(ttf) QCD predictions and subsequent next to leading order (NLO) 
corrections |3] are in good agreement with experimental results, obtained at LHC 
HEIE]. It is clear, on the other hand, that at LHC energies the gluoinc luminosity 
is large, so production of two mesons in two independent partonic interactions 
(so called double parton scattering, DPS) is possible. In papers [3 E] it is shown, 
that cross sections calculated in DPS approximation are comparable with the ones 
obtained in single parton scattering (SPS) approximation. The drawback of DPS 
is poor knowledge of double parton distribution functions. It turns out that while 
for double J/V' production DPS and SPS cross sections are comparable, in the case 
of T(1S')T(1S') hnal state DPS cross section is suppressed by more than order of 
magnitude. 

Question concerning situation in TJ/'ip production arises naturally. At leading 
0{a'^) order in color-singlet approximation this process is forbidden, so the following 
relation holds: 
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( 1 ) 


The reason is that at leading 0 (q:^) order for J/ijjJ/'ilj and TT processes diagrams 
with t-channel quark exchange are allowed, that saturate low invariant mass region 
and lead to rapid decrease of the cross section with the increase of this invariant mass, 
(T ~ 1/s^. At the same order of magnitude the contributions to TJ/'^ production are 
absent, but production via radiative decays of XbXc state is allowed. 

The rest of the paper will be organised as follows. In section [2] the analytical 
expressions for cross sections of partonic reactions gg — )■ Q1Q2 are presented. Cross 
sections of hadronic reactions pp —>■ XbXc -\- X, pp ^ + X are given in section 

[3l The effect of double parton scattering is also discussed in this section. The results 
of our paper are briefly summarised in the last section. 


2 Partonic Cross Sections 

It is well known that production of J/i/jT pair in gg integration in colour-singlet ap¬ 
proximation at leading order is forbidden. Final vector quarkonia particles, however. 
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can be produced via radiative decays of P-wave quarkonia Xb,c, that can be produced 
at LO: 

2 

da{gg J/'ipT + X) = ^ Br {xcj, ^ Br {xbj, ^ Ty) d{gg Xw.XcjJ- (2) 

JbjJc=0 

Typical Feynman diagram of this reaction is shown in FigJU^a), and the corresponding 
matrix element can be written in the form 

M{gg ^ XbJbXcjJ = , Pi — ki)A^^^\k2,pi — fci)e“e 2 + permutations(3) 

In the above expression ei _2 are polarisation vectors of initial gluons with moment 
ki ^2 and we have introduced the notation ^®,(p,k) for the vertex of Q quarkonia 
production in fusion of two (probably virtual) gluons with momenta p, k (see Fig{T](b) 
for the example of the diagram). 
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Figure 1: Typical Feynman diagrams for gg —)■ XbXc process (left hgure) gg —?■ J/'iIjT 
at NLO (right diagram) 


Due to the presented in the amplitude ([3]) propagator of virtual gluons main 
contributions to the cross sections of the considered reaction should come from small 
i kinematic region and in the case of high energy of the partonic reaction this 
gluon can become almost real. In this approximation one can obtain estimates for 
cross sections of scalar and tensor quarkonia production in high energy limit (see also 

H): 


{gg —t XbJbXcJc) 


{2Jc + 1) (2Jb + 1) 


B {XcJc 2g) F {xbj^ —t 2g) 


The ratios of these cross sections are equal to 

d {Xc2Xb2) '■ d {XcOXb2) '■ d {Xc2Xbo) '■ d iXcOXbo) 
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16.3 

14.8 

21.4 

Xh\ 

2.1 

4.6 

3.8 

Xb2 

21.4 

19.6 

29.2 


Table 1: Total cross sections of the partonic reactions gg —)■ XcJcXtJb (iii fb) at s = Sq 


If there is at least one axial meson in the hnal state, the corresponding cross sections 
should be suppressed. 

Presented in papers [El |E1 ES] analysis, however, shows, that experimental data 
are better described with a larger value 

|/?'^(0)|^ = 0.35 ± 0.05 GeV^ 

For bottomonium mesons the same approach gives ini 

|/iM(0)p K 1.0 GeV^ 

These numbers will be used in the following. We will also neglect the hne splitting of 
heavy quarkonia masses, the following values will be used: 


= 3.5 GeV, = 9.9 GeV. 


The scale of the strong coupling constant is chosen to be 


h = 


Xc 


+ M: 


Xb 


Using presented above values it is easy to obtain predictions for partonic cross sections 
of X nieson production in gluon-gluon interaction. Integrated over total phase space 
partonic cross sections of the considered reactions at energy 


s 


So - 2(M^^ -f- 


(5) 


for example, are listed in table [H It is clearly seen that for ratios of scalar and tensor 
mesons production cross sections relations ([1]) are satished with pretty good accuracy. 
Dependence of some of these cross sections on y/S and their distributions over i are 
shown in FigJS] One can see, that in the case of tensor meson production in low t, u 
regions an increase caused by f-channel gluon is observed. For axial meson production, 
on the other hand, in agreement with Landau-Yang theorem such increase is absent. 

In the Introduction it was mentioned already that prohibition of color singlet 
production can be bypassed if one consider NLO corrections (see typical diagram 
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Figure 2: Dependence of cross section of TJ/ip production in partonic subprocess 
(left panel) and its t-distribution at s = sq (right panel). Solid, dotted, and dashed 
lines correspond to XbXc, NLO, and color-octet channels respectively. 




Figure 3: Dependence of partonic cross sections on \/§ (left panel) and i dependence 
at 'hats = So (right panel) for gg —)■ XbXc reactions. On both hgures solid, dashed and 
dotted lines correspond to gg XciXbi, 99 —t XciXb 2 , and gg —)■ Xc 2 Xb 2 subprocesses 
respectively. 
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shown in Fig{T]). In the case of J/ijj pair production this effect was considered in 
paper [13] (see also [H]), and it is simple to extend it to J/V'T case. The cross 
section of the considered reaction can be written in the form 


where 


F{t) 


d&(99 -t TJ/-!/;) 

25vr3a64/2 

'F(i) + F(uy 
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and notations 


( 6 ) 

(7) 


'^b,c 






2k-q 
r = —j—j— 

|q| 


are introduced. One can easily see that after substitution mb = me, et = Cc these 
results agree with those presented in ref.[TH]. The mesonic constants in (I6|) are equal 
to = 0.4 GeV, /x = 0.66 GeV. Using these parameters we have obtained the 
number 


^spf (99 = 21 fb. 

for integrated partonic cross section at s = sq. Dependence of this cross section on 
Vs and its i distribution at s = sq are shown in Figj2]by dotted line. It should be 
noted however, that these results were obtained in the limit s > |tl > SO 

they can be treated only as rough estimates. In our future works we are going to 
remove this limitation and obtain predictions for shown in Fig{T](b) cross sections in 
full kinematical region. 

The other way to bypass prohibition of J/ipT production in gg interaction is to 
consider color octet components of hnal vector quarkonia [12]. Using presented in 
this paper analytical results and 

= 1.3 GeV^ [Oj) = 9.2 GeV^, 

= 3.9 X 10-3 GeV3, (oj) = 0.15 GeV^, 

for NRQGD matrix elements [20l [2111221 [23], it is easy to obtain distribution shown 
by dashed line in Fig[2l According to this hgure y/S dependence of different cross 
sections is quite different, with GO channel giving dominant contribution in region of 
small invariant mass of the hnal pair. 
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Figure 4: Distribution of pp —?■ TJ/'ip+X reaction cross section at LHCb over TJ/'ip- 
pair invariant mass (left figure) and transverse momentum (right figure). Thin solid, 
dashed, dotted, and dash-dotted lines correspond to contributions of XbXo CO, NLO, 
and DPS channels respectively, while thick solid line stands for total cross section 


3 Hadronic Cross Sections 

Let us now consider cross section of hadronic production of J/TT pair in hadronic 
collisions at LHC. Using presented in the previous section results for partonic cross 
sections it is easy to obtain the following values for hadronic ones: 

0'sps(pp-t + ^-t TJ/V'-FX) = 0.2 pb, 

o'sps (pp-t TJ/V'X, NLO) = 1.5 pb, 
c^sps (pp-t TJ/V’,CO) = 11.1 pb. 

It is well known, that at LHC energies double parton scattering (DPS) could also 
give noticeable contribution to cross section of the considered reaction. In this model 
hnal vector quarkonia are produced in two independent partonic subprocesses and 
the cross section of this reaction can be written in surprisingly simple form: 

_ rr tho <^sps(Tf')crsps(<//V’) 

C^DPS ( I d/'?/^) = -. 

^eff 

Using experimental results (Tsps(T) = 0.14 pb, o'spsiJ/'^) = 1-28 pb for SPS cross 
sections and [21] and the value (Te// = 14 mb (HJ [2S1 E], one can obtain the result 

o'DPS (pp-t -1-X) = 12.5 pb. 

It is clear that this simple mechanism gives main contribution to the cross sec¬ 
tion, but one can separate contributions of other channels with the help of various 
kinematical distributions (see. FigsJl|) 
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4 Conclusion 


Presented article is devoted to theoretical analysis of inclusive production of heavy 
quarkonia pair T J/'^ at LHCb. It is clear, that at leading order of perturbation theory 
in colour singlet approximation this process is forbidden, so other mechanisms should 
be studied. In our paper we consider such channels as double parton scattering, next 
to leading order corrections, colour octet model and production of vector quarkonia 
states via radiative decays of Xb,c mesons. For all these processes theoretical predic¬ 
tions of total cross sections under LHCb detector conditions and distributions over 
kinematical variables are presented. 

According to obtained in our article results (see also ra). double parton scattering 
gives main contribution to the cross section of the considered process. From analysis 
of different distributions, however, it is also possible to separate contributions of 
other channel, so one can, for example, determine the value of colour-octet NRQCD 
matrix elements. So, we think that additional theoretical and experimental analysis 
of inclusive TJ/^J pair production at LHC is very interesting and actual task. 

The author would like to thank A.K. Likhoded, S.V. Poslavsky and Dr. Belyaev 
for fruitful discussions. The work was carried out with the hnancial support of FRRC 
and RFBR (grant # 4-02-00096 A) 
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